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I . INTRODUCTION 


The  usual  effect  of  a slight  asymmetry,  caused,  say,  by  an  offset 
center  of  gravity  or  a bent  fin,  is  to  add  a constant  moment  that  rolls 
with  the  missile  to  the  total  aerodynamic  moment  acting  on  the  missile.1 
The  response  to  this  asymmetry  moment  is  a trim  angle  that  also  rolls 
with  the  missile.  The  magnitude  of  this  trim  angle  is  a function  of 
the  spin  rate:  this  magnitude  grows  from  its  value  for  zero  spin,  6, 


TO* 


to  a maximum  resonance  value,  6T_,  reached  when  the  spin  rate  equals 

IK 

the  pitch  rate,  and  then  decays  to  zero  as  the  spin  increases  further. 


We  have  already  studied  the  effect  of  a cubic  static  moment  on 
the  motion  of  a slightly  asymmetric  missile,  using  a general  quasi- 
linear  averaging  technique.2  Except  under  certain  special  spin  condi- 
tions, the  effect  is  to  change  the  frequencies  of  the  transitory  modes 
of  oscillation  and  the  magnitude  of  the  steady-state  trim  angle.  For 
a dynamically  stable  missile,  the  transitory  modes  damp  out  and  only 
the  harmonic  response  to  the  asymmetry  forcing  function  remains. 


Recently,  however,  we  have  shown  that  under  certain  conditions  a 
nonharmonic  steady-state  response  is  possible.3-4  We  call  this  non- 
harmonic response  a generalized  subharmonic  response  because  for  small 
amplitude  motion  it  occurs  at  one- third  the  spin  rate.  Although  the 
theory  of  Reference  3 predicted  two  possible  nonharmonic  motions,  we 
have  been  able  to  generate  only  one  of  these  by  numerical  integration. 
We  conjectured  that  the  other  motion  is  unstable  and  it  is  the  purpose 
of  this  memorandum  report  to  prove  this  conjecture. 


1. 


John  D.  Nicolaides , "On  the  Free  Flight  Motion  of  a Missile  with 
Slight  Configurational  Asyimetries3  " Ballistic  Research  Labora- 
tories Report  No.  858 , AD  26405 s June  1953 ; also  IAS  Preprint  39 5 } 
January  1953. 


2. 


Charles  H.  Murphys  " Nonlinear  Motion  of  a Missile  with  Slight  Con- 
figurational Asymmetries j " J.  Spacecraft  and  Rockets , Vol.  8, 

March  197ls  pp.  259-263 ; also  Ballistic  Research  Laboratories 
Memorandum  Report  No.  20363  AD  870704 s May  1970. 


1 


3.  Charles  H.  Murphy 3 " Generalized  Subhamonic  Response  of  a Missile 
with  Slight  Configurational  Asyrmetriess  " Ballistic  Research  Lab- 
oratories Report  No.  1591s  AD  749787 t June  1972 ; also  AIAA  Paper 
72-97 2 3 September  1972. 


4.  Charles  H.  Murphys  "Subharmonic  Behavior  of  a Slightly  Asymmetric 
Mi8sile3  " AIAA  Journal , Vol.  113  June  1973t  pp.  884-885. 
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II.  NONLINEAR  ANALYSIS 


If  we  limit  the  nonlinearity  under  consideration  to  a cubic  static 
moment  and  neglect  any  linear  Magnus  moment  contribution,  the  transverse 
moment  expansion  assumes  the  form 


1 Ch  = - 1 CM  e1*  + (CO  + c262)i  + CM  £' 
° a 


* CM  u 


lor  zero  spin  and  linear  static  moment  (c2  = 0),  this  moment  expansion 

induces  a trim  angle  of  6j0  . The  ratio  ma  of  the  cubic  part  of  the 

static  moment  to  the  linear  part  for  this  trim  angle  is  a most  impor- 
tant parameter.  It  is  a measure  of  both  the  size  of  the  asymmetry  and 
the  size  of  the  nonlinearity.  By  use  of  this  parameter,  the  nonlinear 
differential  equation  for  the  angular  motion  induced  by  the  aerodynamic 
moment  of  Equation  (1)  can  be  written  as3 


l"  + (H  - iP)l'  - Mq  [1  + ma  (6/6T0)Z]|  = - M()6T0e] 


Equation  (2)  and  various  relations  involving  its  solution  can  be 
considerably  simplified^by  a change  of  the  independent  variable  from 


s to  t,  where  x = (-Mq)  2 s.  If  we  use  dots  to  denote  derivatives  with 


respect  to  x,  Equation  (2)  becomes 


l + (H  - iP )t  + [1  + ma  (6/5t  )2]£  = 6Tne^ 


where  ( ) = ( ) 


The  linear  angular  motion  of  a slightly  asymmetric  missile  (m  = 0) 
is  described  by  the  usual  tricyclic  equation*  a 


l = 


, , i($  + <f>30) 

kj  e + k2  e 2 + k3  e 


Relations  for  the  parameters  ^ ^ , k2  , kg,  $ ) of  Equation 


(4)  can  be  obtained  by  an  averaging  process.2  If  <j>  is  not  near  zero 
or  unity,  these  relations  take  the  form:3 


C*j  - P)  ■ 1 * ma  tSej/5T0>2  i ' 2 


- ■ - i* 


r 


I 


i 


i 


I 


Mr 


f 


_1 


k^kj  = = - [H  \ - 2 ma  k2k3  (sin'r)av][2*1  - P]“  C6) 


k2/k2  = X2  = - [H«j>2  + *2  - ^k3K2‘  isinTjavji.f2 


+ X + m k2k  k:1  (sinY)  1 [2k  - P]"1  (7) 


V = 2<J>1  - 4>2  - 4>  - 4> 


30 


(8) 


k 3 [1  - (1  - I x/!yH2  + ma(2k?  + 2k2  + k3)] 


cos4>  - m k2k  (cos'?) 
30  a i 2 


av 


(9) 


- k3H$  = sin<t>30  + mak2k2(sinH') 


av 


(10) 


where 


6 Aro 

el 


= k2  + 2k|  + 2k2  + 2k2k3(cos4')av 


y4™ 

V / J 


= k2  + 2k2  + 2k2  + kjkgk"1  (cosY) 


av 


In  Reference  3 we  shewed  that  nonharmonic 
can  exist  when  * is  constant.  The  precise  conditions  for  these 

tions  are: 

4/  _ 0 (Ha) 


\1  = x2  = 0 


(lib) 


The  values  of  (kJf  k2,  k3,  $30,  ¥)  associated  with  a steady-state 


solution  will  be  identified  by  bar  superscripts.  Conditions  (11) 
give  a simple  relation  between  kj  and  k2  and  a more  involved  relatio 


between  V and  H : 


11 


, 


XV, 
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(13) 


k2  = 2b  k2 
1 2 

sin?  = (2mak2k3) 


where 

b = - • 1 


Equations  (9,  10,  11a,  12,  13)  constitute  five  equations  for  the  five 
unknown  parameters  (vx,  k2»  k3»  i30.  ?)•  For  the  algebraically  simpler 


case  where  I 
Table  1 . 


/I 

x y 


<<  1,  a concise  form  of  these  equations  is  given  in 


TABLE  1.  CONDITIONS  FOR  STEADY-STATE 

NONHARMONIC  SOLUTION  WHEN 

I /I  <<  1 
x y 

k2  = 2b  k2 
1 2 

(Tl:l) 

sin?  = H4>  [2  (2  + bjmjq^]"1 

(T1 : 2) 

1 + ma[2(b  + l)k|  + 2k2  + 2k2k3cos?]  = [?/(2  + b)]2 

(T1 : 3) 

kg-jl  - ?2  + ma[2(2b  + l)k2  + k2]^-  = cos^30  " 2mab^2  C0S’J' 

(Tl:4) 

sin?3Q  = - H? [k2  + (b/2+b)k2]k31 

(Tl:5) 

<j>x  = £l  + ma[2(b+l)k2  + 2k2  + 2k2k3cos?]J-Js 

(Tl:6) 

?2  = - £l  + ma  [(4b  + l)k2  + 2k2  + 2bk2k3cos'i'iJ-!s 

(T1 : 7) 

b » - b- 

*1 

(T1 : 8) 

12 
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_ iv 
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The  relations  in  this  table  can  be  solved  by  an  iterative  process 
for  a given  m , $ and  H,  provided  that  $ and  ft  satisfy  certain  restric- 


u 

tions.  Equation  (Tl:3)  requires  $ > 3 (that  is,  the  spin  must  exceed 
three  times  the  resonance  spin)  when  m is  positive  and  < 3 when  m 


A 

is  negative.  Equation  (Tl:2)  imposes  an  upper  bound  on  H:  this  bound 

. . A A 


O I ' 

is  shown  in  Figure  1 for  m = ± 0.2  . 

cL 


For  each  H less  than  this  upper  bound,  the  relations  in  Table  1 
yield  two  solutions,  identified  in  Figures  2-5  by  solid  and  dashed 
lines.  Figures  2 and  3 are  plots  of  &2  and  £3,  respectively,  versus 

spin  for  zero  H and  various  m values.  Figures  4 and  5 are  plots  of 


<}>3q  and  Y,  respectively,  versus  H for  m&  = 0.2  and  various  spins, 


The  solutions  identified  by  dashed  lines  in  Figures  2-5  are  the 
ones  that  have  not  been  obtained  by  direct  numerical  integration  and 
the  ones  that  we  want  to  show  are  unstable. 


III.  SOLUTION  FOR  k3  = 0 


As  we  can  see  from  Figure  3,  k3  passes  through  zero  for  one  set  of 

solutions.  The  nature  of  the  solution_is  different  above  and  below  the 
critical  value  of  3>,  say  $ > at  which  k = 0 . Thus  it  is  of  some 

C * 


interest  to  derive  a set  of  equations  for  determining  for  a given 

m . 
a 


From  Equations  (Tl:2)  and  (Tl:5),  we  have 


sin  ^3Q  = - 2mak2|b  k|  + (2+b)k2  J sini 


When  k = 0,  Equations  (14)  and  (Tl:4)  reduce  to 


sin$3Q  = - (2mabk3)sin? 


cos$30  = (2mabk3)cos* 


Thus  we  must  have 


2m  b k?  = ± 1 , k = 0 

3 ^ O 


' .**  ' 


I 


Next,  we  note  that  for  C3  = 0,  Equations  (T1:6-T1:8)  reduce  to 

I2  = 1 + 2(b  + l)m  fc2 
rl  v a 2 


(18) 


<j>2  = 1 + (4b  + l)m  1c2 
2 a 2 


1 + (4b  + l)m  k2 

b2  = 

1 + 2(b  + l)m  k2 

E 4 


(19) 


(20) 


Solving  Equation  (20)  for  m^k2  > we  have 


m k2  = 
a 2 


b2  - 1 
D 


where 

D = D(b)  = 1 + 4b  - 2b2  (1  + b) 
Substituting  Equation  (21)  in  Equation  (17),  we  obtain 

m = 4b2  (b2  - l)3 

^ " D3 


(21) 


(22) 


(23) 


For  zero  k and  any  m , we  can  solve  Equation  (23)  for  b,  then  use 

3 3.  , 

Equations  (18)  and  (21)  to  obtain  the  critical  value  of  <j>  : 


4>c  = (2  + h)4>} 


= (2  + b)  1 + 


[■ 


2 (b+1)  (b2-l) 


r 


(24) 


Critical  values  for  a range  of  m are  listed  in  Table  2 (the  approxima- 

tion  column  is  discussed  in  the  next  section).  Note  that  - 1/27  is  the 
minimum  m for  which  kQ  will  go  to  zero.  As  m increases,  b approaches 

3.  3 3 

bjj  = 1.14,  the  only  positive  root  of  the  cubic  equation 


D = 0 


(25) 


It  can  be  shown  that  whether  or  not  k is  zero,  bM  is  an  upper  bound  on 
b . 3 

14 
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- . . , . ~ : .... 





TABLE  2.  CRITICAL  SPIN  VALUES  (£  - 0) 


m 

a 

b(Eq.  23) 

k2 

1/27 

0.5 

3. 

.03 

0 . 65505 

2.94120 

.01 

0.86737 

3.86298 

.001 

0.95931 

8.04768 

♦c(Eq*  24) 


0. 

0.99684 

1.90778 

2.55635 


0 

1 

00 

3. 

.001 

1.02570 

7.87014 

3.38411 

.01 

1.04567 

3.62960 

3.77834 

.1 

1.07143 

1.67110 

4.51085 

.2 

1.07944 

1.32306 

4.82601 

.5 

1.08962 

0.97180 

5.32761 

1.0 

1 . 09679 

0.76963 

5.78028 

2.0 

• 

1.10333 

• 

0.60965 

6.30441 

• 

• 

00 

• 

1.13973 

0 

00 

0.775 

1.082 

1.890 

2.546 

3. 

3.394 

3.799 

4.549 

4.872 

5.385 

5.848 

6.383 


I,  | | 1- 


IV.  APPROXIMATE  RELATIONS 

As  can  be  seen  from  Figures  2 and  3,  ^ terms  can  usually  be 
neglected  in  comparison  with  k2  terms,  at  least  when  * is  greater  than 
and  not  too  close  to  3.  When  kg  is  neglected,  b is  givenjiy  Equation 
(20) . We  see  from  this  equation  that  for  a given  ma,  as  k2  goes  to 
zero,  b approaches  unity  and  as  k2  increases,  b approaches  the  positive 
root  of  Equation  (25),  namely  1.14.  Thus  a good  first  approxima- 

tion for  b is  1. 

From  Fquations  (8)  and  (11a)  we  have 

7 <j>  T bj>  (28) 

4* i = 2 + b * 92  ~ ' 2 + b 

Hence  for  b = 1,  our  first  approximation  of  the  frequencies  is: 

♦ * - i i C29) 


Then  from  Equations  (18-19,  29),  we  have 

*1  s 1 + 4raa^  s ir 


♦1  = 1 + 5ma^2  * 9“ 


Equations  (30)  and  (31)  yield  two  possible  approximations: 

m C2  r <b2  - 9 ~ ?. 

ma  2 " 36  ’ 45 

The  average  gives  an  excellent  first  approximation  for  k^* 

» H i ^ 1 

a 2 40 

Substituting  this  result  in  Equations  (30-31)  we  have  the  improved 
frequency  approximations: 


: 'I’fcit'i  i ‘ •,iin~te  i ivni1!  .I  ... 


16 


from  which  we  can 


obtain  an  improved  approximation  for  b. 


~r5(f  - ni  *.r, 

[4  (<J. 2 + 1)J  L 4(*  + 15J 


a ! . - 9 ■ gi!  - 1- 

s'ci’2  + it  + iy 


Simple  approximations  for  k3>  4>30  and  V can  now  be  obtained. 
First,  we  neglect  k2  in  Equations  (Tl:4)  and  (14)  so  that  they  become, 


respectively, 


2mabk|  cosi  - cos<|>3q 
E3  5 4»2  - 1 - 2(2b+l)mak^ 


sin*  = - 2m  bk^  siny 
30  a 2 


= - |m  \~h  f (^sinY 
a 1 1 


where 


= 2b(mak|)|mak2r 


with  b and  m k2  given  by  Equations  (35)  and  (32).  Next,  we  substitute 
the  approximation  (36)  for  k3  in  Equation  (Tl:2)  to  obtain,  with  the 

help  of  Equation  (37) : 


' 

' . >■ 


ft*  [t2  - 1 - 2 (2b  + l)mak2] 


2 (2  + b)m  k 
v a 2 


^2mabk|  cos'?  - cost^  sin* 


= - sin$„_  cos*  - cost  sin? 
y30  30 


= - sin(t30  + *) 


where 


H = - Ha  sin(t30  + 


2 (2  + b)maE2 

tp2  - 1 - 2 (2b  + l)mak2j 


m f (<j>) 


2(2  + b)  lma^|l^ > Q 

= t[t2  ‘ 1 - 2(2b  * l)mak2J 


and  where  b and  mak2  are  again  given  by  Equations  (35)  and  (32) . Thus 

the  elaborate  damping  curves  of  Figures  4 and  5 would  be  transformed 
into  a family  of  near  sine  waves  if  H were  plotted  against  4>30  + V • 

Note  that  H.  has  the  sign  of  m . Hence  approximation  (39) 
predicts  that  when  H is  at  its  maximum  value  (£  |HA|),  then 


d>  + 
*30 


T =/’/2,  »a  < 

^3tt/2,  ma  > 


(42) 


We  will  test  the  adequacy  of  this  prediction  in  the  next  section. 


The  process  of  approximation  cap  be  summarized  as  follows.  For  a 
given  m , $ and  H: 


1.  compute  b and  m k2  , Equations  (35)  and  (32); 


2.  compute  fJf  f2,  HA  and  * + ^ , Equations  (38),  (41),  (40) 
and  (39) ; 


3.  from  f and  5 + ^3Q  , compute  and  f , Equation  (37); 
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4.  compute  k , Equation  (36). 


An  approximation  for  the  critical  spin  $c  can  also  be  easily 

obtained.  Using  Equation  (32)  and  setting  b = 1,  we  can  write 
Equation  (17)  as: 


_4_ 

m 


or 


$2  = 9 


(43) 


Critical  spin  values  computed  from  Equation  (43)  are  given  in  the  final 
column  of  Table  2.  Note  that  although  Equation  (43)  can  be  applied  to 


ma  > - 4(9/40) 3 = - .046,  the  approximation  is  poor  below  ma  = - .03. 


1 


MM 


V.  EXACT  SOLUTION 

In  Figure  6,  exact  solutions  of  Equations  (Tl:l-5)  are  given  in 

the  form  of  5 versus  ¥ plots  for  ra  = ± 0.2  and  various  values  of  $ . 
30  r a 

Note  the  two  families  of  curves  for  m = 0.2  . For  spins  between 

cl 

3 and  the  critical  value  of  4.826,  ¥ grows  from  0 to  i while  $3Q  grows 

from  n to  (<j>  ) < 3tt/2  and  then  decreases  back  to  tt  . The  value  of 

30  MAX 

(ijj  ) can  be  estimated  by  setting  ¥ = it/2  in  Equation  (37).  For 
30  MAX 

spins  greater  than  the  critical  value,  grows  from  ir  to  2tt  while 
¥ grows  from  0 to  (¥)MAX  < tt/2  and  then  decreases  back  to  0 . The 
value  of  can  be  estimated  by  setting  $30  = 3ir/2  in  Equation  (37) 

For  m = - 0.2  , there  is  only  one  family  since  no  critical  spin 

3- 

value  exists.  The  equations  of  Table  1 lose  their  validity  as  we 
approach  resonance  (<fi  = 1);  hence,  curves  for  spin  values  less  than  2 
are  not  shown. 

When  H = 0 (and  i<3  / 0) , Equations  (Tl:2)  and  (Tl:5)  reduce  to 

sin¥  = sin^30  = 0 (4^) 

For  any  value  of  and  ma,  the  particular  values  of  (¥,  3Q)  for  the 

two  possible  zero-damping  solutions  are  given  in  Table  3 . 


TABLE  3.  ZERO  DAMPING  POINTS 


neg 

neg 

(-"> 

2tt) 

(-IT,  IT) 

neg 

pos 

( o. 

it) 

(-IT,  IT) 

pos 

neg 

( o. 

tt) 

( IT,  IT) 

pos 

pos 

( o. 

it) 

( 0,  2 it) 

%**4- 


k 


' ' ?-*  ■ 


■ 


" TFT 


The  labeling  of  the  points  as  stable  or  unstable  is  at  this  point  in 
the  report  only  a supposition.  In  Figure  6,  of  course,  only  the  last 
three  pairs  in  Table  3 occur. 


As  we  move  from  (0,  it)  along  any  curve  4*  = constant  in  Figure  6, 
H increases  until  it  reaches  its  maximum  value  approximately  at  the 
intersection  of  that  curve  with  the  line  L^: 


LA  : *30  * 


. r»/2, »,  < o 

|3tt/2,  m&  > 0 


(45) 


That  is  (as  we  noted  in  the  previous  section),  line  approximates 
the  maximum  damping  locus 


^ : H ~ 


(46) 


which  is  our  conjectured  stability/instability  boundary.  To  see  just 
how  well  La  approximates  LH<  we  computed  the  exact  values  of 


U H (*30  ♦ f)fi 


hmax 


(47) 


for  a variety  of  <j>  and  m values.  Figure  7 plots  the  difference 


U - 3ir/2  versus  $ for  m =0.2  and  0.5  . It  was  found  that  for 


positive  m , this  difference  is  very  nearly  proportional  to  rtn 

cit 


a 


U - 3tt/2  » vST  F($)  > 0 

ci 


(48) 


Note  from  Figure  7 that  the  approximation  5 LH  worsens  as  we 

approach  $ = 3;  this  is  because  out:  assumptions  in  deriving  approxima- 
tion (39) (in  particular,  that  k3  could  be  neglected  in  several  equa- 
tions) are  inadequate  near  ♦ » 3 . However,  for  * greater  than,  say, 
five,  L.  is  an  excellent  approximation  to  our  conjectured  boundary. 

/V 
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VI . STABILITY  ANALYSIS 


In  order  to  determine  the  stability  of  a nonharmonic  steady-state 
motion  identified  by  the  parameters  (fcj,  &2»  &3»  *30»  we  must  con' 

sider  the  behavior  of  motions  near  such  a limit  motion.  Therefore,  we 
must  derive  the  differential  equations  for  small  perturbation  function 
ti.(t)  defined  by  the  equations 


ic.  + n. 

(49) 

l l 

k + n 

(50) 

2 2 

it  + n 

(51) 

3 3 

$ + n. 

(52) 

y + n 

5 

(53) 

If  Equations  (49-53)  are  substituted  in  Equation  (5),  the 
frequency  *.  (j  . 1,  2)  can  be  expressed  as  the  sun  of  the  frequency 

J.  for  the  limit  motion  (t,.  t,,  .*)»»>»  linear  fu"Cti°" 

J - • . . 


<t> 

of  the  perturbation  variables: 


4.  * 4.  + 6.  j = 1 *2 
3 J 


(54) 


where  the  e.  = ry  y ng)  are  given  in  Table  4 for  « 1 

Differential  equations  for  \ and  n£  can  now  be  obtained  by 
substituting  Equations  (49-53)  in  Equations  (6-8);  these  too  are  given 
in  Table  4 . 


22 


TABLE  4.  PERTURBATION  EQUATIONS  FOR  I /I  « 1 


«,  ■ cvvLv,  * C2\  * N cos*)n2 

+ (2k ^ + k2COS?)n3  - (k2k3sin^nl 


(T4:l) 


3 2 3 


= 0 n /♦  ) (2b)*  (2k  + k cos^ni 


+ (k  - bk  cos4’)n 

v 2 3 2 


(T4 : 2) 


+ (2k  + bk  cos^n  - (bfc  k sin^n, 

' a n O 9 4 


\ - - (v2*}  « ei  + ki 


- 2m  f(k  n + k n )sin?  + H n cos V 
a(_  32  23  235 


\ - - (v2*|fi 


(T4:3) 


e + e 
2 2 


♦ ^[tao*  k3nj  - bk3n2  * bk2n3jsin'» 

* a*.bi2t3V01*  I 

err  A . 4 


(T4:4) 


fi3  . - 2% 


- 2maj(2k3  + k^cos^k^  + (2k3  + bk^cos’VJk^n^j 

- H ^ ♦ jV  - 1 - »a  (2k2  * 2k2  ♦ 3k23)]n3 


ri,  = - 2m 


+ 2k  Yn  - (sin 5 )n  + (2m  bk^sin^n 

3 4 30  4 a ^ 3 

_ f A 

- 2m  k2sin»  (k^  + bk2n2)  - $(2n3  + HV 

- H kan..  - (cos^30)n4  - (Zn^bkJcos'On,. 


(T4:5) 


\ = 2ei  - e2  ~ \ 


k - i ^ r , a ^ 


; T/  "'V  ' , 

• *V  V . 


IMi 


.... . 


...  ...  .. — ^w... 


The  differential  equations  for  n3  and  n^,  however,  must  be  obtained 

from  Equation  (3).  In  our  previous  use  of  Equation  (3),  the  phase  and 
amplitude  of  the  one-  and  two-modes  were  allowed  to  be  functions  of  the 
independent  variable  in  the  assumed  quasi-linear  solution.  Equation 
(4),  but  k and  <j>3Q  were  required  to  be  constant. 


In  order  to  derive  the  proper  equations  for  ng  and  r»4  , our  assumed 


quasi-linear  solution.  Equation  (4),  must  be  enlarged  to  include  these 
variables:  r , . i f a * 7 + n 


? = 


+ k2e 


+ (k3  + n3)e 


♦ A ♦ n ) 

*30  h 


Equation  (55)  is  now  substituted  in  Equation  (3);  the  result  is  divided 
by  exp  i(<}>  + <T3Q  + n4)  and  averaged  over  a distance  that  is  large  com- 
pared to  the  various  wavelengths  involved.2  The  real  and  imaginary 
parts  of  the  resulting  equations  are  the  differential  equations  in  n3 
and  n given  in  Table  4 . 


The  usual  procedure  for  stability  analyses  is  to  assume  solutions 
to  the  perturbation  equations  of  the  form 


n3  = V 6 


1,  2...N 


where  N,  the  number  of  parameters,  is  in  our  case  five.  Equation  (56) 
is  then  substituted  in  the  N perturbation  equations;  the  result  is  a 
system  of  N homogeneous  linear  equations  in  the  hjo’s,  with  coefficients 


involving  X . Let  DET  be  the  N-th  order  determinant  of  this  system. 
If  M is  the  sum  of  the  orders  of  the  N perturbation  equations  (M  = 7 
for  the  equations  of  Table  4) , the  equation  DET  = 0 yields  an  M-th 
order  polynomial  equation  in  X,  the  characteristic  equation  of  the 

system: 


VM  * v"’*  * • • • * Vi1  * ‘ 0 


where  the  a^s  are  assumed  real.  Let  r^  r2  •**  rM  be  the  r00tS 
of  this  equation.  Then  the  system  represented  by  Equation  (57)  is 


1.  stable  if  the  real  parts  of  the  M roots  are  all  negative: 


R(rk)  < 0,  k * 1,  2,  . . . M 


, V ...  \ 

}•.  * - ♦ j*  . 

: . • ■ : 


2.  critically  stable  if  (57)  has  one  or  more  non-repeated  pure 
imaginary  roots  (including"the  root  0 i) , but  the  remaining  roots  have 
negative  real  parts; 

3.  unstable  if  (57)  has  either 

a.  one  or  more  roots  with  positive  real  part 


b.  repeated  pure  imaginary  roots. 

Note  from  2 and  3 above  that  if  the  characteristic  equation  (57)  has 
one  zero  root,  the  system  will  be  either  critically  stable  or  un- 
stable,  but  if  (57)  has  more  than  one  zero  root,  the  system  can  only 
be  unstable. 

It  is  not  necessary  to  evaluate  the  roots  to  test  for  stability, 
most  stability  tests  are  made  directly  on  the  coefficients  . For 

zero  damping  (see  the  next  section),  fairly  elementary  tests  suffice; 
for  nonzero  damping  (see  the  final  section),  we  will  require  the  more 
sophisticated  Routh-Hurwitz  stability  criterion. 

VII.  STABILITY  FOR  ZERO  DAMPING 

Numerical  integrations  of  Equation  (3)  led  us  to  suspect  that 
the  zero-damping  point  (0,  it)  in  Figure  6 and  Table  3 represents  un- 
stable solutions.  To  prove  this,  we  set  H = 0 (and  hence  sin<|>30 
sin'F  = 0)  in  the  perturbation  equations  of  Table  4 and  make 
the  simplifying  assumptions  listed  in  Table  5.  The  resulting  appro 
imate  perturbation  equations  (T5:l-7)  are  also  listed  in  Table  . 

Substituting  the  assumed  solution  (56)  in  Equations  (T5:l-7),  we 
obtain  the  characteristic  equation: 

X7  + a2X5  + a^X3  + agX  = 0 C59) 

When  a k,  a4  and  a6  have  been  simplified  (by  the  same  assumptions  made 
in  Table  5 and  by  the  use  of  Equation  (Tl:4)  with  b - 1,  kg  = 0),  we 
obtain: 

a2  = 2(1  + 4>2  + 6m&k2) 

a4  = (1  - *2  + + 72Cma^2)2 

- H7ma2k3  [(1  + 34>2) k 3 + cos^30Jcosi/<P 
a6  = 9ma2R3|l3(4»2  - 1)  - 60mak2  ] cos^cosV/l2 


8 /•  'i  " ;*■ 


* 


where 


TABLE  5.  APPROXIMATE  PERTURBATION  EQUATIONS 

FOR  I /I  <<  1 AND  fi  = 0 

x y 


ASSUMPTIONS: 

b = 1 

ic3  « max  £1,  k2) 


k.e.  <<  h.  j B 1»  2 

J 3 3 

l = - L - */3 
1 2 


EQUATIONS: 

ej  = (3mak2/*)[/2  nj 


+ 2 n2  + (cos9)n3 


e2  = ‘ L2/2  ni  + n2  + (cos^T,s 


n,  = ^2  (3mak2k3cos'i1/*)n 


n = (3ra  k2k  cos¥/*)n_ 

2 a 2 3 5 


^ = - 2m  k^cosYCv^  n.  + O 

^ a 2 * ^ 


(T5 : 1) 
(T5:2) 

(T5:3) 
(T5 : 4) 


+ (4>2  - 1 - 6ma^^n3  + 2V\ 

(T5: 5) 

ri4  = - 2$n3  - (cosj^n,  - (2mak*cos*)n5 

(T5:6) 

ng  = (3mak2/i)  |^4v^  O 5 ^ + (3cosY)n3j 

- n (T5:7) 
4 

cos*  = ± 1 
30 

cosY  * ± 1 

’ 1 W" 1 


1M 


m 


or,  using 


Approximation  (32)  and  assuming  4>  > 2.5, 


a2  5 (23/10)*2 


a,  = (17/20)  r 
4 


= (207/80)  (d>2  - 9)mak2cosi3ocos4’ 


(60a) 


(60b) 


(60c) 


we  obtain  the  sixth  order  equation 

(X2)3  * »2P>2)2  + a*cx2)  * \ ■ 0 t61) 

. , . , -I  of  this  cubic  in  X2  yield  the  remaining 

The  three  roots  (Sj,  s2>  s3J  o 

six  roots  of  Equation  (59)  . 

A = 0,  ± , ± ^ ± *^3 

■ • i c+aHn l i tv  s s and  s,  must  be  real,  negative  and 

For  critical  stability,  Sj,  s2  a 3 

unequal : 

1. . H woulf  LveUarpos?five^ealUre 

imaginaries;  otherwise  +vs^  si 

component; 

2.  unequal  because  repeated  pure  imaginary  roots  indicate 
instability. 

' tiaet  c sn  and  Sq  be  real  and  unequal  is  satis- 

The  requirement  that  sj,  s2  ana  53  nepative: 

fied  if  the  discriminant  of  Equation  (61)  is  negative. 

q3  + < 0 (62) 

where  3a4  - a^ 


9a  a - 27a  - 2a3 

D - 94  £ - 

R “ 54 


■ 

, V 


5V1’' 


where 


R = 


37513  t, 
216000 


0.1744>6 


Hence  Inequality  (62)  is  always  satisfied: 

q3  + r2  S - 0.0116<5>12  < 0 

and  the  three  roots  are  rea]  and  unequal. 

To  determine  when  the  three  roots  are  negative,  we  use  the  Descartes 
Rule  of  Signs:  the  number  of  positive  real  roots  of  the  cubic  (61)  is 
equal  to  (or  two  less  than)  the  number  of  sign  changes  in  the  sequence 
1,  a2,  a4,  a0  . Since  a2  and  a4  are  positive,  there  will  be  no  sign 

changes  and  hence  no  positive  real  roots  if  and  only  if 

a > 0 (63) 

6 

Since  we  have  already  shown  that  all  three  roots  are  real.  Inequality 
(63)  is  the  condition  for  three  negative  roots. 

To  summarize:  for  zero  damping,  critical  stability  requires  that 

Sj,  s2  and  s3  be  real,  negative  and  unequal;  the  three  roots  are  always 

real  and  unequal;  they  are  negative  if  and  only  if  condition  (63)  holds. 

Now  consider  our  expression  for  a , Equation  (60c).  Since  (d>2  - 9)m 
is  always  positive  and  since  cos<t>30  = ±1»  cost  = ±1  for  zero  damping, 
condition  (63)  reduces  to 

cos^QCOsi1  = 1 (64) 

That  is,  for  zero  damping,  our  system  is  critically  stable  when  (64) 
holds  and  is  unstable  when  the  product  of  the  cosines  is  -1. 

This  proves  our  conjecture  on  the  instability  of  the  point  (0,  11 } 
and  justifies  the  labeling  of  the  other  points  in  Table  3 as  stable  or 
unstable.  We  can  now  turn  to  the  stability  analysis  for  nonzero 
damping  with  some  confidence  in  our  approach. 


VIII.  STABILITY  FOR  NONZERO  DAMPING 


For  nonzero  damping,  the  assumed  solutions  (56)  were  substituted 
in  the  perturbation  equations  of  Table  4.  None  of  the  assumptions 
and  approximations  that  simplified  our  zero-damping  analysis  were 
applied  to  the  non zero -damping  equations.  As  a result,  the  coeffi- 
cients of  the  seventh-order  characteristic  equation: 


X7  + a X6  + a X5  + a,X4  + a X3  + a,.X2  + agX  + a?  = 0 (65) 

1 2 3 


were  truly  horrendous  functions  of  m&,  and  H (none  of  them  vanishing 

everywhere  as  did  a^  a3  and  a5  for  zero  damping).  The  opportunities 

for  elementary  algebraic  blunders  were  many  and  we  seemed  to  miss  few 
of  these  opportunities.  However,  by  a series  of  independent  checks 
and  double-checks  we  obtained  what  we  managed  to  convince  ourselves 
was  an  error-free  set  of  expressions  for  the  a^s.  To  test  these  co- 
efficients for  stability,  we  used  the  Routh-Hurwitz  criterion  given  in 
Table  6. 


TABLE  6.  ROUTH-HURWITZ  STABILITY  CRITERION 


The  real  parts  of  the  M roots  of  the  equation 


VM  + aTxM  1 + •••  + ®M-1  X + ®M  ■ ° 
will  all  be  negative  if  and  only  if 


(T6: 1) 


T.  > 0 , i = 0,  1,  2,. ..M 

l 


(T6:2) 


where 


T = a 
0 o 


(T6:3) 


Tl  - ai 


(T6:4) 


and  T n i 2,  is  the  n-th  order  determinant  of  the  upper  left 
n 


nxn  elements  of  the  array 


(T6:5) 


. . . , ^ 


■MR 


’ - — ■ 


Our  complete  procedure  was  as  follows.  For  a given  m^  and  <j>  : 

A _ — . 

1.  We  assigned  a value  to  any  one  of  the  trio  H,  <j>3Q>  Y and 

solved  (by  an  iterative  process)  the  equations  of  Table  1 for  the  re- 
maining parameters.  That  is,  we  obtained  a solution  point 

5=  Ck,.  k2,  k3,  $30,  *) 

• ys  o 

for  a specified  m&  and  <|>  and  a specified  or  determined  H.  (The  option 
to  evaluate  results  at  a specified  $3Q  or  V » rather  than  at  a specified 
H,  proved  useful:  in  certain  regions  of  the  4*,  $3Q  plane,  the  con- 

vergence of  the  iterative  process  was  improved  by  fixing  one  or  the 
other  of  the  two  angles.) 

2.  At  point  A we  evaluated  the  coefficients  a^  of  the  character- 
istic equation  (65)  formed  from  the  perturbation  equations  of  Table  4. 

3.  Finally,  we  evaluated  the  Routh-Hurwitz  determinants  of 
Table  6 to  test  point  A for  stability. 

The  above  process  was  incorporated  into  a single  computer  program 
and  applied  to  hundreds  of  input  values  so  that  we  could  map  out  the 
stability  regions  in  the  solution  plane  and  define  the  stability/in- 
stability boundary  numerically. 

• A A , 

For  a given  m^  and  <|>  and  for  any  H less  than  two  solutions 

were  obtained:  one  stable  and  one  unstable.  In  every  case,  the 

boundary  occurred  (within  the  accuracy  of  the  computations)  at  a point 
where  H = . This  was  precisely  our  conjecture. 


30 


• 0.2 


r**-. 


CM  iqQ 


Figure  5.  H versus  ? (m  - 0.2,  various  *) 
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the  locus  of  solution  points  at  which  4>30  + Y = 3ir/2; 
approximately  L^ 


/V  A 

the  locus  of  solution  points  at  which  H = Hj^;  the 
conjectured  stability/ instability  boundary 

reference  length 


c <5 
2 TO 


CM 

a 

4>  =0 

CM 

a 

<S2  = 0 

order  of  the  characteristic  equation  (57) ; 
sum  of  the  orders  of  the  perturbation  equations 


PS&3 
21  C0 


number  of  parameters 
I 


-p-  t the  gyroscopic  spin 

y 


X,  Y,  Z components  of  the  missile's  angular  velocity 
components  of  the  cubic  discriminant.  Equation  (62) 
k-th  root  of  the  characteristic  equation  (57) 


reference  area 

arclength  along  trajectory,  in  calibers 
roots  of  the  cubic  equation  (61) 


i-th  Routh-Hurwitz  determinant.  Table  6 
the  value  of  (cj>  + Y)  at  H = H^x  for  £ 


the  value  of  (^^  + ^)  at  H = H^x  for  a given  <f> 
and  m 

a 


i ' 'i  - It*  ‘'J 
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LIST  OF  SYMBOLS  (continued) 

Y,  Z components  of  the  velocity 
magnitude  of  the  velocity 

nonrolling  Cartesian  coordinate  axes.  The  X-axis 

lies° along  the  missile's  axes  of  ^7t^V-axis 
Y and  Z axes  are  so  constrained  that  (1) 
is  initially  in  the  horizontal  plane  and  (11)  the 
angular  velocity  of  the  coordinate  system  has  a 
zero  X- component. 

|jr|,  the  sine  of  the  total  angle  of  attack 
trim  angle  at  resonance 
trim  angle  at  zero  spin 


£1  ’ S 


perturbation  in  4>  , 4>2 

value  of  at  t = 0,  Equation  (56) 


V V V V n5  Perturbation  in  V V V *30’  ’ 


damping  factor  in  the  assumed  solution.  Equation 
(30)  , to  the  perturbation  equations 


(q  + ir)E 

V 


v + iw 
V 

air  density 

(-Mo)*4  s 
roll  angle 


spin  = roll  rate  = p 

critical  spin;  for  a given  ma,  the  4>  value  at 


+ ...  ... ^ if..  , . 

l*~"  r “ 1 - ' nw«i»tg|  i*  ■,  "f/r'jT 
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LIST  OF  SYMBOLS  (continued) 


1-  and  2-mode  phase  angles,  Equation  (4) 


^30 


initial  value  of  the  trim  mode  phase  angle,  Equation 

(4) 

2<p  - <f>  - <f>  - > Equation  (8) 

1 2 u « 


Superscripts 

c y 

n 

C) 


d(  )/ds 

d(  )/dT 

c ) 


c ) 


value  associated  with  some  nonharmonic  steady-state 
motion 


n 


value  in  the  nonrolling  coordinate  system 


Subscripts 

( ) average  over  a distance  that  is  large  compared  with 

the  wavelengths  involved 

( ) maximum  value  for  a given  <j>  and  m 
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